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We consider the phenomenological description of the two pion invariant mass 
spectrum in the V — > Vtttt decays. We study the parametrization of the 
amplitude involving both S and D wave contributions. From a fit to the two 
pion decays of the T(nS) and ^(nS) we determine the /o(600) mass and 
width to be m fo = 528 ± 32 MeV and T fo = 413 ± 45 MeV. The mass and 
width values we report correspond to the real and imaginary part of the S 
matrix pole respectively. 



The experimental identification of low mass scalar resonances is a long 
standing puzzle whose origin can be traced back to some of the following 
characteristics: a large decay width, possible mixing with multiquark or glue 
balls, overlap of resonances and the opening of channels, responsible for the 
appearance of cusp effects [1]. In particular the /o(600) has a long history, 
it has been included in some issues of the Review of Particle Properties but 
it has also been excluded for long periods by the Particle Data Group. Re- 
cently, experimental evidence from different corners of particle physics has 
accumulated confirming the existence of the /o(600) resonance. There have 
been attempts to interpret the low lying scalars as multiquark states [2, 3] or 
KK bound states [4, 5]. Also, models exist starting directly from chiral sym- 
metry [6, 7, 8] or else unitarized models where the scalar nonet arises [9, 10]. 
However the understanding of actual processes involving scalar mesons for 
descriptions starting from first principles are not available, effective theo- 
ries are not well suited to deal with these scalars, the use of sum rules is 
perhaps the best approach to the problemfll], however no predictions for 
all of the scalars have been advanced in that framework. Furthermore, the 
phenomenological description of broad resonance faces severe problems, the 
determination of the physical parameters -mass and width- is a non trivial 
problem that requires study of the non-resonant background dependence. 

For small invariant mass of the pion pair, data is dominated by the phase 
shift [12, 13]. Studies have been carried in this kinematical region using 
chiral symmetry and Roy equations -solid theoretical tools- establishing thus 
a firm result to be considered by other analysis [14]. Inclusion of the di-pion 
low invariant mass favor a light and very broad / (600) (m ps 470 ± 30MeV^ 
and T/2 ps 295 ± 20 MeV). Below 1 GeV, information on the tttt phase 
shift is obtained from ttN scattering, PP annihilation at rest and central 
production, these data allow the existence of a broad (r ps 500MeV) scalar 
meson resonance [15, 16]. Decay of pseudoscalar charmed mesons are also 
a source of valuable data involving scalar mesons, although the / (600) has 
only been reported by the E-791 collaboration after the Dalitz plot analysis of 
D — > 7r + 7r~7r + [17]. Another important source of information are the vector 
meson decays. The processes considered involve either the scalars themselves 
or pion pairs together with photons or vector mesons. Among these we can 
mention (p, u, $) -> P Pj [18, 19, 20, 21, 22], J/V -> P P V , -> P P, 
and T(nS) — > T(mS)P P [23], where P stands for a pseudo scalar (n or K) 
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and V for a vector meson (p,ou). Experimental evidence for the contribution 
of scalar resonances to some of these processes has been reported, and a recent 
analysis of the data for T(nS) — > T(mS)P P concludes the contribution 
of the / (600) with a large uncertainty in the width (m = 526I37, V = 
SQl^^M eV) [26]. It should also be noticed that in the experimental data 
for ip(2S) — > 7r + 7r~ J/^> reported by the BES collaboration [25], no evidence 
for / (600) contribution is foreseen. 

In this paper we concentrate on the decay of heavy quark vector meson 
resonances (T, and where pair of pions are produced with invariant mass 
ranging below the 1 GeV region. The kind of processes we are interested in 
has been considered by a number of authors using techniques as diverse as 
pure chiral symmetry [27], non-relativistic theory assuming the existence of 
an Adler zero [28], the color field multipole expansion in the non- relativistic 
limit [29, 30], effective Lagrangian based on chiral symmetry and the heavy 
quark expansion [31, 32, 33], and also a purely phenomenological description 
based on a Breit-Wigner parametrization [26]. The framework so developed 
is used then to describe the two pion invariant mass spectrum and, in some 
cases, also the angular distribution. The latter is important since existing 
experimental data can discriminate among the proposed models. Further- 
more, as far as we know, the complete parametrization of the amplitude has 
not been discussed and some confusion has arisen concerning the "S" and 
"D" wave contributions to the amplitude [33]. 

Our purpose is to perform an analysis as general as possible, including 
both S and D waves in the (PP) system. We show that Lorentz covariance 
fixes the parametrization of the amplitude in terms of four form factors. Us- 
ing the two gluon mechanism implied by the OZI rule, the identification of 
the spin zero (di-pion S wave) and spin two (di-pion D wave) are unambigu- 
ously identified. For simplicity we restrict our analysis to two form factors, 
ao(^7T7r) and a2(m 7Tn ) which are associated to S and D waves respectively. 
We then invoke the pole approach to parametrize ao, a<i- We propose a Breit- 
Wigner plus a background for a and a soft background for a 2 . In this way we 
claim that crossed channel as well as higher scalar resonance contributions 
are taken into account. Thus, we expect that the amplitude associated to 
crossed channel and higher resonance contribution behave as a soft function 
in the 500-800 MeV range, where the / (600) is expected to lie. Our strategy 
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is to perform a joint fit to data from different processes, involving 165 points 
and 33 parameters. 

The set of data points we consider include the tt — tt invariant mass 
distribution of the T(3S) -> T(1S) + tttt [23], T(3S) -> T(2S) + tttt [23], 
T(2S) -> T(l 1 S) + 7T7r [23, 24], ^(25) -> J/V' + tttt [25] decays. The following 
characteristics are worth-mentioning: i) we are considering flavor conserv- 
ing processes, and all of them are expected to proceed through two gluons. 
This point will be relevant when discussing the parametrization of the form 
factors, ii) the smallness of the phase space available for the processes un- 
der consideration < y/s < 0.9 GeV). Note that the expected central 
value and the large width of the / (600) would imply non negligible reso- 
nance effects in these processes and in) the invariant mass distribution for 
s — * s th -where s t h stands for the threshold value of the di-pion invariant 
mass- show a peculiar behavior to be contrasted with the typical (s — 2m 2 ) 
expected in processes involving soft pions. Compare for example in Fig(l) 
the threshold behavior of the yfs distribution for T(2S') — > T(IS') + tttt or 

1p(2S) -> J/lp + TTTT With T(3S) -> T(1S) + TTTT. 

In order to understand the nature of the problem we face, remark that the 
more recent data [24, 25] have been analyzed in terms of the scale anomaly. 
Indeed, these processes can be fitted without problem, which brings the ques- 
tion if the full set of data we consider can be explained using the same for- 
malism. Our results show that this is not the case and that inclusion of the 
/o(600) improves our understanding of the data. Thus, as far as we can see, 
any attempt to provide a successful description of the flavor conserving two 
pion decays of the T and tp families should consider the full set of data, since 
the foreseen physical mechanisms (scale anomaly, scalar resonance exchange) 
are operative in all cases under consideration. 



1 Parametrization of the amplitude 

We are interested in the decay V'(p',rf) — > V(p, ^)vr(pi)7r(p2) where the 
letters in parenthesis stand for the four-momenta and polarization of the 
corresponding particles. We introduce q = pi — p 2 and Q = p\ +p 2 - Q 2 = s, 
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and p' 2 = m' 2 ,p 2 = m 2 ,pi 2 = P2 2 = m 2 . In order to obtain the general 
parametrization is convenient to consider the amplitude for the exchange 
of arbitrary spin zero and spin two meson like objects (although we do not 
consider the actual exchange of any particle). Let us first consider the S wave 
contribution. The amplitude for the V — > V+ scalar can be written as: 

M = 77 W 

Covariance allow us to write t^ u in terms of g^ v and the independent (p',p) 
four-momenta. Imposing p' • rj' — 0, p • rj — and using p'H^ = 0, which 
follows from the fact that V'(p') is produced through a virtual photon in an 
e + e~ machine, we obtain: 

Mo = a Q 77 • rj 1 

V vp J 

Already at this point we encounter differences with the parametrizations 
used in the literature, where only the rj ■ rj' term is considered. Although 
sizable effects are not produced by the extra term it is important to work 
with the proper Lorentz invariant amplitude. In particular, differences could 
become relevant when polarization measurements are involved. 1 For the D 
wave contribution we write: 

M 2 = A^U^B^ (2) 

where A^ u describes the V'(p') — > V(p) + D(Q), D standing for a spin two 
meson like object, B pa describes the D — > nn amplitude and II^ pcr is the 
spin two projector: 

5 

Rfivpa = h^{\)h p(T {\). (3) 
A=l 

The polarization tensor has the following properties: 

iNote that q 2 J2 po i l M o| 2 = 2 + ?f> whcreas J2 po i(v " »?') = 2 + S> with m and E thc 
outgoing vector meson mass and energy respectively. 
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M*) = M*)> = g^K^X) = o (4) 

Using these relations together with the projector property of Yi^ upa one finds: 



with 



P„u = 9pu - —qT- ( 6 ) 

On the other hand, by Lorentz covariance B pa oc q p q a - It is convenient to 
introduce: 



IV = U^ pa B^ = q p q v - l -P pu q 2 . (7) 

Furthermore, using Lorentz covariance and imposing again the conditions 
p' ■ r)' = 0, p ■ r) = and p'^A^ = 0, due to the coupling of V' with a 
virtual photon, we find the general structure for A pu , and through Eq. (2) 
the di-pion D wave amplitude (Note that Q^Yi pv = Q u n pu = 0, therefore 

P nv = ^iv). 

M 2 = b -j^jW- p)) n ^ + c (v'^v ■ p') - p' ^ v ' p ^' p) ) v lv ^v 



" 2 ^,')-^)«, (8) 



m' 2 — m 2 \ p • p 



Where we introduced the m' 2 — m 2 factor to work with a dimensionless a 2 . 

In order to obtain the decay rate, we carry out the integration over the pi, pi 
Lorentz invariant phase space in the two pion center of mass reference frame, 
i.e. Q — 0, q = 0. We obtain: 
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r 1 1 [\^\\st 12 rf3 P rf3 ^l ^2 , n4,4/^ n 



- 1 /'vim 2 rf3p n *"V 

The S-D wave interference vanishes upon integration. Indeed covariance 
imply 



but using Q^IT^ = g^H^ = 0, it follows that a = b = 0. 
For the .D wave we proceed along the same lines. 



1 1 y f d 3 p (A IJ - v A* pa Mi 
2m' 3 fjJ (2vr)3(2p ) 1 ' ' 



where 



7 P oZ (2vr) (2pi ) (2tt) (2p 20 ) 

= xIT 

x is determined using the reference frame where Q = 0, go = 0. For example 

TT3333 = 2 r = f E(n33n 33 ) - ^ - - , (2tt) V(g - Pl - Pa ), 

3 J po ; (2tt) (2pio) (2tt) (2p 20 ) 

in this way we obtain 
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60tt 



_ 4m| 

Q 2 ) 



We observe that using the q = 0, Q = reference frame, it is easy to 
show that 11^ ~ Y 2 (9,(f>), i.e. it is associated to the di-pion D wave. This 
together with the fact that a 2 , b and c can only depend upon s = (pi + p 2 ) 2 , 
allow us to conclude that M. 2 describes the di-pion spin 2 wave. In the 
following we consider the particular case where only one Lorentz invariant 
amplitude is included in A^ u . To this end we set b = c = in Eq(8). This 
is the simplest way to consistently introduce D wave effects in the invariant 
amplitude. Within this approximation we obtain for the di-pion invariant 
mass distribution (s = {p 1 + p 2 ) 2 = Q 2 ). 



dY {m'p)^fs 
~d\ft ~ 3(4m'7r) 3 



pol 



(y-p')W -p) , 
p ■ p' 



{SW + DW) 



(9) 



where p stands for the three momentum (p = |p|) and 



SW=\a (s)\ 2 (l 



Am 



(10) 



DW 



02(f) 

n' 2 ~m' 2 

180 



:i- 



Am 2 



s — (m! + m )) — 4m' m 



(11) 



with 



2m p= (s — m — m) — 4m m 



(12) 



pol 



(r] ■ p')(r}' ■ p) 



p.p> 



= 2 + 



4m 2 m 



/2 



(m 2 + m' 2 — s) 2 



(13) 
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We have now a general expression describing the decay V — > Vtitt, which 
involves two invariant amplitudes, associated to the S and D wave respec- 
tively. We assume, when y/s < 0.9 GeV, the S wave is composed of the 
/o(600) and a non resonant background. This is an approximation since 
other resonances could contribute to the amplitude in the kinematical re- 
gion considered. However, it is reasonable to expect that the contribution 
of the / (980) (ry ( 980 ) ~ 100 MeV) and higher states decaying in two pi- 
ons behave softly in the neighborhood of the / (600), even if the latter is a 
broad resonance. Crossed channel contributions are treated in a similar way, 
i. e. considered as soft functions of the di-pion invariant mass. Therefore we 
parametrize the form factor associated to the S wave in the following way. 



II In 



itlq is introduced for dimensional reasons and is fixed to itlq = 0.5 GeV. For 
D(s) we used two different expressions: 



D(s) = s-m 2 + U{s) 

D(s) = s — m 2 p + im p T ' p (15) 

where n(s) stands for the /o(600) self energy which involves loop of kaons and 
pions. Note that Eq(15) defines the mass m p and width T p of the resonance 
in terms of the real and imaginary part of the pole of the S matrix. If a pole 
exist, it should be independent of the process where it is observed. However 
neither the residue nor the background have to be the same for different 
processes. For this reason we include the index % which is associated to the 
physical decay under consideration. In the kinematical region of interest the 
D channel is non-resonant, thus we can approximate it by a soft background: 



a 



(0 

2 I i _ « (; 



(16) 



It should be clear that our approach is a phenomenological, and that we have 
not attempted to explicitly incorporate the scale anomaly. In fact, it can be 
argued that such a contribution is included within the background. 
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2 The fit 



We consider the following decays: T(3S) -> T(1S) + tttt [23], T(3S) -> 
T(2S) + tttt [23], T(2S) -> T(IS') + tttt [23, 24], ^(25) -> J/^ + tttt [25], 
resulting in a total of 165 points. All the data points but the last have been 
gotten from plots since no listings of the data are available. It is important 
to mention that previous fits [26] of the ip(2S) — > J/ip + tttt included only a 
subset of the experimental data. 

Following the analysis in [24, 25] we first attempted a fit in terms of the 
scale anomaly [29, 30]. The fit produces a Xd.o.f. > 2. This is not an 
unexpected result, the amplitude associated to the scale anomaly is derived 
under the assumption that the pions are soft, and this is not the case for 
all the data under consideration. In fact, according to [26] and previous 
analysis [1], the / (600) is expected to contribute to these processes. This is 
the motivation to try a fit in terms of the pole approach, as discussed in the 
previous section. 

The fit using the pole approach is based on Eqs(9,10, 11,14,15,16). The 
pole can be parametrized in terms of the full propagator, including the self- 
energy II(s), which involves loops of kaons and pions. In such result 
of the fit, we get a mass and the gjKK an d gjmr coupling constants, and from 
these the pole position is determined. Unfortunately, the fit turns out to be 
insensitive to the values of the coupling constants, i.e. different values of the 
g's lead to essentially the same fit, but correspond to pole positions far from 
each other. For this reason our analysis is restricted to D(s) given by Eq(15) 
where the mass (m p )and width (r p )are obtained directly from the fit. 

The fit involves 33 parameters. We consider four processes (i — 1 to 4 in 
Eqs(14, 16)) and for each of these we require seven parameters (b and / which 
are complex, and a, c, g). Three normalization factors are also parameters of 
the fit -because some of the reported data refers to number of events, not to 
a differential decay rate- and finally the mass and width of the resonance. 
The fit leads to a pole located at m p = 528 ± 32 MeV and T p = 413 ± 45 
MeV together with the parameters reported in Table(l), and a^L/ = 112. 
The result of the joint fit are shown as the solid lines in Fig(l). We do 
not report the result of the fit for S wave alone, however we should mention 
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that including D wave effects improves the total x 2 but the Xdof remains 
unchanged. Note in this respect that the fit parameters are obtained from the 
di-pion invariant mass spectrum, the angular distributions are not involved. 
Once the fit parameters are fixed, numerical integration over the di-pion 
invariant mass is performed and we obtain -up to normalization again- the 
following angular distributions (in GeV units) and the curves shown in Fig(3): 



T{2S) -> T(l 1 S)7T7r[24] : — = 0.07(0.0674(3cos 2 # - l) 2 + 1.09) (17) 



d(cos9) 



JT"1 

T{3S) -> T(l 1 S)7T7r[23] : — — = 2.2(0.0009(3cos 2 # - l) 2 + 0.511) (18) 

CL I COSu J 



V(2S) -> J/^tvtv[25] : — = 6.5(5.49(3cos 2 # - l) 2 + 359) (19) 

a(cos9) 

Besides the pole position and the quality measure through the Xd.oj.i ^ 
is instructive to analyze the different contributions to the decay rate. These 
are shown in Fig(2) where the dots corresponds to the pole contribution, 
the dashed line to the background and resulting fit is represented by the 
continuous line. From these plots one can see that the pole and background 
contributions must interfere destructively in all cases, except the T(3S f ) — > 
T(1S) transition for which the kinematical region includes the pole position 
and both, constructive (below the pole-mass) and destructive (above the 
pole) interferences occur leading thus to the structure around 650 MeV. In 
table(2) we quote the contributions to the x 2 from the different processes. 



3 Summary 

In this paper we consider the phenomenological description of the follow- 
ing decays: T(3S) -> T(1S) + tv + tv, T(3S) -> T(2S) + tv + tv, T(2S) -> 
T(1S) + tt + 7r, ip(2S) — > J/ip + 7i + TV. As far as we can see, it is not 
possible to obtain a good quality fit in terms of the scale anomaly alone. 
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Using general arguments, we derived an expression for the invariant ampli- 
tude describing flavor conserving processes of the type V — > V + n + n, 
including both S and D waves, which involves two invariant amplitudes. We 
parametrized the S wave form factor with a pole plus a soft background and 
the D wave form factor by pure soft background. Fitting the data yields a 
pole in m p = 528 ± 32MeV and T p = 413 ± AhMeV with xlo.f. = 112 - 
We remark that a strong interference among the pole and the background is 
required to fit the data. Thus, our analysis seems to indicate that physics in 
S wave pion-pion interaction below 800 MeV is governed by a subtle inter- 
play between _/o(600) meson contributions and a big background, difficult to 
understand in terms of conventional physics and which could be associated 
to the scale anomaly. 
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FIGURE CAPTION. 



Figure 1. Data points used in the analysis [23, 24, 25] and the resulting fit 
(solid line). In the horizontal axes we plot m nn = y/s whereas the vertical 
refers to the differential decay rate J^, or number of events, as shown in the 
plots and discussed in the quoted references. For the same reaction, open 
circles and solid triangles refer to data obtained from exclusive and inclusive 
processes respectively. 

Figure 2. The curve resulting from the fit to the data (solid line), and for 
comparison the pole (dashed line) and background (dotted line) contribu- 
tions. Note that for the T(3S) — > T(2S) + tttt the fit-curve is close to the 
axes, which implies a strong destructive interference among the pole and the 
background contributions. 

Figure 3. Angular distribution as obtained from Eqs(17,18,19) and compared 
to data from Ref([24, 23, 25]). Open circles and solid triangles refer to data 
obtained from exclusive and inclusive processes respectively. 



TABLE CAPTION. 



Table 1. Parameters resulting from the fit. The normalization factors refer to: 
N a data from Ref([24]), N 2 to the T(3S) -> T(15)+7r°7r° and N 3 to T(3S) -> 
T(2S) + 7T 7r°. The parameters are defined by Eqs(14, 16). Processes are 
labeled 1: T(2S) -> T(15)7T7t, 2: T(3S) -> T(1S)tvtt 3: T(3S) -> T(2S)nn 
and 4: V(2S) -> J/^nn. 

Table 2. Xd.o.f. f° r eac ' 1 separate process as obtained from the fit to the 
di-pion invariant mass distribution. 
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T(2s)->T(1s) TT + 7T(Ref. 23) 



T(2s)>T(1s) ttV (Ref. 24) 
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Figure 3. 
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N 2 
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N 3 


0.62 




8.2xl0 3 
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5.3xl0 3 + i5.1xl0 3 
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-0.14 
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92 
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Table 1. 



Process 


■Ndata 


xio. f. 


T(2s) -> T(1S)7T7T 


48 


0.74 


T(3S) -> T(1S)7T7T 


40 


1.3 


T(3S) -> T(2S)7T7T 


32 


1.1 


■0(2s) — ► J/ipnn 


45 


1.5 



Table 2. 
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